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HESS J1825— 137 was detected with a significance of 8.1 a in the Galactic Plane survey conducted with the 
H.E.S.S. instrument in 2004. Both HESS J1825-137 and the X-ray pulsar wind nebula G18.0-0.7 (associated 
with the Vela-like pulsar PSRB1823-13) are offset south of the pulsar, which may be the result of the SNR 
expanding into an inhomogeneous medium. The TeV size (~ 35 pc, for a distance of 4 kpc) is ~ 6 times larger 
than the X-ray size, which may be the result of propagation effects as a result of the longer lifetime of TeV 
emitting electrons, compared to the relatively short lifetime of keV synchrotron emitting electrons. The TeV 
photon spectral index of ~ 2.4 can also be related to the extended PWN X-ray synchrotron photon index of 
~ 2.3, if this spectrum is dominated by synchrotron cooling. The anomalously large size of the pulsar wind 
nebula can be explained if the pulsar was born with a relatively large initial spindown power and braking index 
n ~ 2, provided that the SNR expanded into the hot ISM with relatively low density (~ 0.003 cm -3 ). 

1. Introduction 

PSRB1823-13 is a 101 ms evolved pulsar with a spin-down age of T = 2.1 x 10 4 years (for an assumed 
braking index of n — 3) |4 1 and in these properties very similar to the Vela pulsar. It is located at a distance of 
d = 3.9 ± 0.4 kpc 10. High resolution XMM-Newton observations of the pulsar region confirmed a previous 
ROSAT detection of an asymmetric diffuse nebula, which was hence given the name G 18.0-0.7 For the 
compact core with extent i?cN ~ 30" (CN: compact nebula) immediately surrounding the pulsar, a photon 
index of Tqn = 1-6^Q 2 was measured with an unabsorbed luminosity of Lqn ~ 9<i| x 10 32 ergs -1 in the 
0.5 to 10 keV range for a distance of Ad^ kpc. The corresponding pulsar wind shock radius is R s < 15" = 
0.3(^4 pc. The compact core is embedded in a region of extended diffuse emission which is clearly one-sided, 
revealing a structure south of the pulsar, with an extension of _Ren ~ 5', (EN: extended nebula) whereas the 
~ 4' east-west extension is symmetric around the north-south axis. The spectrum of this extended component 
is softer with a photon index of Ten ~ 2.3, with an unabsorbed luminosity of Len = 3^4 x 10 33 erg s _1 for 
the 0.5 to 10 keV interval. No associated supernova remnant (SNR) has been identified yet. 

2. A possible association of HESS J1825-137, G18.0-0.7 and 3EG J1826-1302 

At 7-ray energies, PSRB1823-13 was proposed to power the close-by unidentified EGRET source 3EG J 1 826- 
1302 |6]. The region around PSR B 1823-1 3 was observed as part of the survey of the Galactic plane with the 
H.E.S.S. instrument |1|. In this survey, a source of very high-energy (VHE) 7-rays (HESS J 1825-1 37) 11' 
south of the pulsar was discovered with a significance of 8.1 a. We note, that the new VHE 7-ray source is 
located within the 95% positional confidence level of the EGRET source 3EG J1826-1302 and could therefore 
be related to this as of yet unidentified object. More details of the H.E.S.S. observations, analyses and merits 
of the possible association of HESS J1825-137 with G 18.0-0.7 was discussed by 0. 

Figure 1 shows a north-south declination slice of G18.0-0.7 in X-rays and HESS J1825-137 in 7-rays. Both 
images show that the emission peaks at the position of PSR B 1823- 13, but with extended emission shifted 
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towards the south. The prominent X-ray synchrotron peak at the pulsar position is as a result of the compact 
X-ray nebula, where the magnetic field strength should be largest. Whereas the extent of the X-ray emission is 
only ~ 5 arcmin, the TeV emission is seen out to ~ 0.5°. Even though the position of the unpulsed EGRET 
source 3EG J1826— 1302 is consistent with that of HESS J1825-137 |2|, the poor angular resolution of the 
EGRET instrument precludes a similar morphological study. 
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Figure 1. (a) left: north-south declination slice of G 18.0-0.7 in X-rays 1 8 1 showing the bright compact nebular core at the 
pulsar position and the ~ 5' extended X-ray nebula extending to the south of PSRB 1823-1 3 (indicated by a vertical line), 
(b) right: similar slice for HESS J1825-137, showing a similar extension south of the pulsar, but with a much larger size of 
~ 0.5°. 



3. The one sided nature of the PWN and the anomalous size of HESS J1825-137 



The one-sided nature of this PWN as seen in TeV was already suggested by |8|, for the X-rays, based on the 
hydrodynamical simulations of |3| for Vela X and earlier studies referenced by |8|: It was assumed that the 
density of the interstellar medium (ISM) around the progenitor star was inhomogeneous along the north-south 
direction, with the density /pressure towards the northern direction significantly larger relative to the south. The 
reverse shock from the northern direction should then have crashed relatively early into the PWN, pushing the 
latter towards the south, as observed. 

The remaining problem with the association of HESS J1825-137 with G18.0-0.7/PSRB1823-13 is the anoma- 
lously large size of the TeV source: the 0.5° diameter of the source translates to a radius of Rpwn ~ 17^4 
pc, and in the absence of radiative cooling in a SNR shell, one would expect a SNR forward shock radius of 
Rsnr ~ 4:Rpwn — 70d,4 pc |93. One possible key towards solving this problem would be to consider a 
relatively large initial spindown power (Lq) and pulsar braking index less than the canonical value of n = 3, 
which would increase the age T of the SNR. Integrating the spindown power Epsr = J dt lilCl over the 
increased age would also maximise the PdV work done by the PWN against the SNR. The Sedov-Taylor SNR 
size would also increase with increasing age (scaling as T 2 / 5 ), giving a radius of 
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Here E§n ~ 10 51 ergs is the SN explosion energy and N ~ 0.003 cm -3 is the density of the hot phase of 
the ISM |7|. A pulsar braking index of n = 2 ("LMC pulsar" type) will increase the age of PSRB 1823-1 3 to 
~ 42 kyr if we assume a relatively large Lq (i.e. an initial spin period much less than 0.1 s). This will imply a 
size consistent with the predicted size of ~ 70 pc. 

Figure 2(a) shows the evolution of EpgR with time for different values of n and L : it is also clear that the 
total PWN energy is maximised for minimum n and maximal value of L , but it is not yet clear how to relate 
this time dependent spindown power to the evolution of the PWN in a quantitave way. 
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Figure 2. (a) left: integrated spindown power for PSRB 1823-13 assuming different braking indices n and initial spindown 
power Lq values as indicated. The arrow indicates the PdV work required to displace the observed PWN volume against 
a pressure of 5 eV/cm 3 . The vertical position of the arrow would then scale linearly with pressure, (b) right: the spectral 
energy distribution of the PWN of PSRB1823-13: HESS J1825-137 (triangles), G 18.0-0.7 (compact and extended nebula 
indicated by two lines) and 3EG J1826— 1302 (circles). For details see 0. 



4. The spectral energy distribution of the PWN of PSR B1823-13 

Aharonian et al. [2| derived the spectrum of HESS J1825-137, which can be represented by a power law in 
energy, with photon index 2.40±0.09 stat ±0.2 sys and a flux above 230 GeV of (3.4±0.2 stat ±0.8 sas ) x 10" 11 
cm _2 s _1 (corresponding to 12% of the Crab flux above that energy). This spectrum is shown in Figure 2 (b), 
together with the X-ray spectra of G18.0-0.7 GO and 3EG J 1826-1 302 0. It is clear that the TeV spectrum 
of HESS J1825-137 does not extrapolate towards GeV energies, and a spectral flattening below ~ 10 GeV is 
required. This spectral break may be the result of a cooling effect, which is also seen in X-rays: The compact 
nebular X-ray photon spectral index of ~ 1.6 represents the uncooled component, whereas the photon index 
of ~ 2.3 measured from the extended X-ray nebula, represents the synchrotron cooled component (as a result 
of propagation effects). The TeV photon index of ~ 2.4 is then consistent with inverse Compton scattering 
of these synchrotron cooled electrons on the CMBR, provided that this cooled spectrum also extends into the 
EUV/soft X-ray band, where electrons contribute directly to the observed TeV emission for a field strength of 
~ 10/iG or less. 

The expected intrinsic break in the injected spectrum, represented by the so-called 'T-factor" of the wind at the 
pulsar wind shock of the compact nebula may also contribute to the observed flattening: in this case the number 
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spectrum of electrons with energy E e > Tm e c 2 is expected to follow an accelerated power law spectrum of 
the form E~ 2 - 2 , resulting in the observed (uncooled) compact nebular spectrum with a photon index of 1.6 (as 
discussed above). Electrons with energy E e < Tm e c 2 is then expected to have a much harder spectrum. This 
component is however unconstrained by radio or optical observations. 

We also note that the energy flux in TeV is larger than the synchrotron energy flux of the extended X-ray neb- 
ula. This dominance of the TeV flux holds even if we extrapolate the spectrum of the extended X-ray nebula to 
the EUV/soft X-ray band. The evaluation of this soft synchrotron component is important, since electrons con- 
tributing to this soft unseen component, also scatter CMBR photons into the H.E.S.S. band for field strengths 
of ~ 10/iG or less: Note that the XMM X-ray spectrum is only measured within 5' from the pulsar, whereas 
electrons contributing to this soft band should have a much longer lifetime. They can therefore propagate to 
larger distances (as seen in TeV), resulting in a population of accumulated electrons over the lifetime of the 
pulsar, whereas the X-rays seen in Figure 2 (b) represent synchrotron emission from freshly injected electrons. 
This accumulation effect should increase the intensity of the soft component accordingly. Thus, the soft un- 
seen synchrotron component should have a similar size as seen in TeV. A test for this hypothesis would be to 
measure the size of G 18.0-0.7 in the softest channels of XMM over an area similar to the TeV size. A more 
detailed study is however required to model the intensity of this unseen soft synchrotron component. 

5. Appendix: lifetimes for X-ray and TeV emitting electrons 

The respective synchrotron lifetimes of electrons in a transverse magnetic field of strength B = 10~ 5 £?_5 G, 
radiating synchrotron photons with mean energy i?keV (in units of keV) and inverse Compton scattered CMBR 
photons to a mean 7-ray energy of i/rw ( m units of TeV), are given by 

r x = 1.2 kyrBZ 3 5 /2 E-^ 2 and r 7 = 4.8 ky^I^E'^ 2 . 
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